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1.  INTRODUCTION.  The  applicability  of  the  unconstrained  adjoint  vari- 
ational statement  in  solving  nonconservative  stability  problems  has  been 
shown  in  a series  of  articles  [1-4].  The  problems  involved  the  stability 
of  beams  or  columns  subject  to  concentrated  or  distributed  tangential  loads. 
The  problems  are  solved  by  finding  the  variational  statement  associated  with 
the  differential  equation  and  they  are  rendered  unconstrained  by  incorporating 
the  geometric  boundary  conditions  into  the  variational  statement  through  the 
use  of  Lagrange  multipliers.  With  the  variational  statement  now  available, 
the  problem  is  discretized  and  solved  using  finite  elements.  Various  types 
of  external  forces  and  geometric  boundary  conditions  can  be  handled  using 
the  above  techniques.  It  is  the  purpose  of  this  paper  to  incorporate  into 
the  above-mentioned  formulation  the  effect  of  support  damping  and  to  examine 
its  effect  on  the  solution.  The  specific  problem  chosen  was  that  of  a 
rotating  cantilever  beam  of  constant  cross  section. 

This  problem  was  chosen  for  its  application  to  a simplified  helicopter 
blade  and  although  no  nonconservative  forces  are  considered  the  solution 
technique  outlined  above  is  applicable. 

The  effects  of  support  damping  on  the  vibration  response  of  beams  has 
been  investigated  by  others.  Fu  and  Mentel  [5]  and  Mentel  [6]  considered 
support  damping  due  to  viscoelastic  layers  applied  to  the  ends  of  a beam 
in  its  supports.  The  effect  of  translational  (axial)  damping  was  found  to 
be  of  the  same  order  in  terms  of  energy  dissipation  at  the  supports  as  that 
of  material  damping.  Material  damping  effects  were  found  to  stiffen  the 
beam  which  increased  both,  the  resonance  frequency  and  the  energy  dissipa- 
tion at  the  supports.  They  also  found  that  rotational  motion  dominates  the 
damping  properties  if  all  parameters  are  suitably  optimized.  Ruzicka  [7] 
presented  an  evaluation  of  the  resonance  characteristics  of  undirect ional 
vibration  isolation  systems  including  directly  coupled  (Kelvin/Voight)  and 
elastically  coupled  (Zener)  damping  elements.  His  results  were  mostly  for 
the  Zener  model  and  he  found  that  resonant  frequencies  of  vibration  isola- 
tion systems  with  viscous’ damping  may  increase  or  decrease  with  an  increase 
in  the  viscous  damping  coefficients  depending  on  the  stiffnesses  in  the 
system.  MacBain  and  Genin  treated  support  flexibility  in  a series  of 
papers.  Support  and  material  damping  was  introduced  in  [8].  The  support 
is  viewed  as  a complex  rotational  support  stiffness  based  on  bounds  for 
the  elastic  modulus  found  in  their  earlier  papers.  They  find  that  when  the 
support  damping  constant  is  an  optimum,  the  support  loss  factor  is  also  an 
optimum,  and  system  loss  factor  reaches  a maximum  value.  This  same  value 
of  the  support  loss  factor  is  also  that  which  critically  dampens  the  system 
in  free  vibration. 

The  results  presented  here  show  the  effects  of  support  damping  on  the 
flexural  frequencies  of  vibration  of  the  rotating  beam  with  both  deflection 
and  rotation  flexibility  at  the  support. 
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2.  PROBLEM  STATEMENT.  The  geometry  of  the  problem  is  shown  in  Figure  1, 
The  beam  has  a constant  cross  section  of  area  A,  density  p,  Young's  modulus, 

E,  and  moment  of  inertion,  I.  The  beam  rotates  about  an  axis  fixed  at  one 
end  of  the  beam  and  is  flexibly  supported  at  that  end  by  a deflection  spring, 
kj,  and  a rotation  spring,  k2*  Viscous  dashpots,  c^  and  C2,  are  assumed  in 
parallel  to  the  deflection  and  rotation  springs,  respectively.  The  beam 
rotates  at  constant  angular  velocity,  ft.  S(0)  represents  support  reaction. 


The  differential  equation  governing  the  motion  is  given  by  [9] 

U""  - [ (£2  - x2)U']'+^-u=  0 

2EI  1 El 

and  the  boundary  conditions  are 

at  x = 0 

C2  k2 

u" (0)  - — u'(0)  - — u'(0)  = 0 
El  El 


(1) 


at  x = £ 


C-*  k, 

u'"  (0)  + — u(0)  + u(0) 

El  El 


u" (£)  = 0 
u"'  (£)  = 0 


^51  u<  (0)  = 0 

El 


(2) 


(3) 


The  differential  equation  and  boundary  conditions  are  rewritten  using 
dimensionless  variables  and  parameters  defined  by  the  following: 


u = — 


_ x n _ ftzpA£‘t  ; _ r HI  t1/2» 
T * x “ — * *•  “ l ■;  „ttJ 


2EI 


pA£ 


Cl  = 


Cj£ 


(EIpA) 


kl  " 


,1/2 

(EIpA)1/2£ 

(4) 

kx£3 

k2£ 

- 

El 

2 EI 

Time  is  removed  by  assuming  displacements  to  have  the  form 

u(x,t)  = u(x)e^ 

Then  the  differential  equation  becomes  (dropping  the  bar  symbol) 

u""  - Q[ (1  - x2)u']’  ♦ X2u  = 0 


(5) 


(6) 


and  the  boundary  conditions 


r u" (0)  - (Xc2  + k2)u'  (0)  = 0 
1 u"'  (0)  + (XCl  + k^uCO)  - Qu'(0)  = 0 


(7) 


The  eigenvalues,  X, 


u"(l)  = 0 
u'"  (1)  = 0 
will  be  complex, 

X = Xp  + iX  j 


The  frequencies  will  be  given  by 


u)  ~ X j 


[-51-] 1/2 

pAr 


(8) 


(9) 


(10) 


and  for  this  problem,  Xp  < 0,  i.e.,  the  real  component  of  the  eigenvalue 
is  negative  and  no  instabilities  should  exist. 


3.  VARIATIONAL  STATEMENT.  To  find  the  form  of  the  variational 
statement,  the  differential  equation  is  multiplied  by  an  arbitrary  variation 
of  the  adjoint  field  variable,  <5v(x),  and  integrated  over  the  beam  length. 
Integration  by  parts  indicates  the  form  of  the  variational  statement  and 
the  natural  boundary  conditions.  The  geometric  boundary  conditions  are 
attached  with  the  values  of  the  springs  and  dashpots  playing  the  role  of 
Lagrange  multipliers.  The  variational  statement  is  finally  given  by 


where 


6J  = 0 


j = /V 
0 


+ Q [ 1 - x2]u'v'  + X2uv]dx  + 


(11) 


+ (Xc1  + k1)u(0)v(0)  + (Xc2  + k2)u' (O)v' (0)  Cl  2) 

Performing  the  variation  of  J with  respect  to  u and  v,  one  can  arrive 
at  the  original  boundary  value  problem  as  well  as  its  adjoint.  In  this  case 
the  two  problems  are  identical. 

4.  FINITE  ELEMENTS.  To  solve  the  problem  using  finite  element  tech- 
niques, the  beam  must  be  divided  into  segments  and  the  nodes  defined.  The 
value  for  the  unknown  variable  within  each  element  must  then  be  expressed 
in  terms  of  the  nodal  values  of  the  function  through  the  use  of  interpo- 
lating shape  functions.  A global  expression,  or  matrix  is  then  formed  and 
the  eigenvalues  found. 


3 


f 


The  procedure  begins  by  taking  the  variation  of  Equation  (12)  and 
allowing  the  variations  in  the  problem  variable,  <5u(x),  to  be  zero, 

/1[u"«v"  + Q(1  - x2)u'6v'  + X2u6v]dx  + 

0 

(Xc2  + k1)u(0)6v(0)  + (Xc2  + k2)u' (0)6v' (0)  = 0 (13) 

The  beam  is  divided  into  L elements,  letting 


£ = L {x  - i = 1,2,3. ..L 

L 


(14) 


be  the  running  coordinate  in  each  element.  Substituting  Eq.  (14)  into  Eq. 
(13): 


l /1[L3u(l)M6v(l)"  + r (L2  - [£  + (i  - l)]2)u(l-),6v(i)' 


i=l  0 

I 2 


+ £ u(i)6v(l)]dS  + (Xcj  + k1)u(1)(0)6v(1)(0)  + (Xc2  + k2) 

L2u(1)'(0)<5v(1)'(0)  = 0 (15) 

In  order  that  the  displacements  and  their  derivatives  within  an  element 
be  expressed  in  terms  of  their  nodal  values,  the  coordinate  vectors 
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[j(i)T  = {u^  U2^  U^J 

and  . (16) 

v(i)T  . (»{»  v2(i)  v«  v™) 

are  introduced,  , U2^  represent  the  displacement  and  slope  at  the 

left  end  of  the  ith  element  and  and  represent  deflection  and 

slope  at  the  right  eijd.  A similar  interpretation  is  applied  to  the  adjoint 
coordinate  vector  v(1).  The  transform  is  indicated  by  T. 

Hermitian  polynomials  are  used  to  relate  the  displacements  within  an 
element  to  its  nodal  values,  hence,  the  following  shape  function  is  assumed, 

iT(S)  ■ (1  - 3£2  ♦ 2£3  £ - 2£2  + £3  3£2  - 2^  -£2  ♦ £3}  (17) 
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So  that 


u(l)(£)  = aT(C)U(l) 


v(i)(5)  = iT(5)vW 

Substituting  Eq.  (18)  into  Eq.  (15) 


X u(i)T{l3c  + [QL  - ? (i  - 1)2]B  - ? E - 


i=l 


where 


2 (i  - 1)  i 6 + 


+ (Xcj  + k1)U(1)TH6V(1)  + L2(Xc2  + k2)U(1-)TF6V(1 


A = /1a(C)aT(C)dC  , E = /Va* (^)aT’ (5)d5 


B = f^'CO iT'(5)d5  , F = a'  (0)a,T(0) 


C = J1 a"(C)aT"(()d( 


D =/^I'(QiT'(C)dC  , H = S(0)iT(0) 


Regrouping  of  (19), 


ZG^VpW  + XRfl)  + S(l)}6V^l)  = 0 

where 

P(l)  = A/L  i = 1,2, ...L 

R(1)  = + CjH  + c2FL2  i = 1 

R(l)  =0  i = 2,3, ...L 

S(1)  = L3C  + QLB  - Q/L  E + kjH  + l^FL2  i 

S(l)  = L3C  + QL[1  - f2(i  - 1 ) 2 ] B - £ E - 2(i  - 1)QD/L 


Using  certain  continuity  conditions  between  the  element  nodal  values 


,(i)  (i-l  (i)  (i-ll 

L = U_  V,  = V J 

13  1 ’ 

I(i)  = U(i_1)  v(i)  = v^'1) 
2 4 2 4 


(25) 


One  can  write 

U(T){X2[P]  + X[R]  + [S] }6V  = 0 (26) 


where  now 


-m 

■ <u'l> 

,(!) 

(D 

u(1> 

U(2) 

uf2>...u(L) 

UL 

U.  } 

1 

2 

3 

4 

3 

4 3 

4 

vT 

- (V(1) 

V(D 

v(1) 

V(D 

v(2) 

(2)  (L) 
V1"  ;...V 

vtL)} 

1 

2 

3 

4 

3 

4 3 

4 

[P],  [R].  [S] 

are  N x N 

matrices  (N 

= 2L 

+ 2). 

Since  6V  is 

arbitrary,  the 

eigenvalue  problem  reduces  to 

U(T){X2[P]  + X [ R]  +[S]«  0 (27) 

for  the  eigenvalues  of  the  problem. 

An  existing  subroutine  was  used  to  find  the  eigenvalues  which  required 
the  standard  eigenvalue  problem  form 

{[A]  + X[I] }U  = 0 (28) 

The  equation 

{X2[A]  + X[B]  + [C] >U  = 0 (29) 

can  be  reduced  to  Eq.  (28)  by  defining 

W = XU 

This  leads  to  the  matrix  equation 

_[0] | [IJ_ 

-[A]_1[C]  | - [A] -1 [B] 

which  is  in  the  required  format.  The  drawback  here  is  that  the  order  of 
the  matrix  has  been  doubled.  Equation  (31),  however,  is  the  form  used  for 
computing  the  eigenvalues. 
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5.  RESULTS.  Figures  2 and  3 show  the  effects  for  zero  damping  at 
the  support.  Figure  2 shows  the  effect  of  the  rotation  spring  (k2)  only 
on  the  frequency  with  load  as  a parameter.  The  deflection  spring  is 
assumed  to  be  infinitely  stiff.  For  Q = 0,  the  beam  is  only  vibrating 
and  is  not  rotating.  One  can  see  that  a stiffening  effect  occurs,  i.e., 
the  vibrating  frequencies  increase  with  an  increase  in  the  rotation  spring 
constant.  The  frequencies  rapidly  approach  those  for  a fully  clamped  vibrating 
and  rotating  beam.  These  results  also  fall  within  the  bounds  computed  by 
Boyce,  DiPrima  and  Handelman  [10].  In  Figure  3,  the  rotation  spring  is 
assumed  to  be  infinitely  stiff  and  the  effect  of  varying  the  deflection 
spring  is  shown  for  different  loads.  Again,  in  general,  there  is  a stiffening 
effect  as  the  deflection  spring  value  increases.  For  very  small  values  of 
the  deflection  spring,  the  first  vibrating  frequency  decreases  slightly  for 
increased  loads,  although  only  Q = 0 and  200  are  shown.  As  kj  is  increased, 
there  are  cross  over  points  after  which  higher  loads  do  imply  higher 
frequencies. 

Figure  4 shows  the  effect  of  rotation  damping  at  the  support  of  a beam 
having  rotation  flexibility  at  the  support.  The  deflection  spring  is 
assumed  infinitely  stiff  and  the  deflection  dashpot  is  zero.  The  figure  is 
for  a specific  value  of  the  rotation  spring,  k2  = 1,  and  shows  the  first 
two  eigenvalues  for  each  of  two  loads,  Q = 0 and  Q = 100.  A stiffening 
effect  is  found  for  increasing  damping  for  Q ^ 0.  For  Q = 0,  there  is  very 
slight  decrease  for  very  small  damping  values.  For  smaller  rotation  spring 
constants,  and  zero  load  frequencies  decrease  with  increased  damping  as 
shown  in  Figure  4 by  the  portion  of  the  results  for  k2  = . 1.  These  results 
are  better  shown  in  Figure  5.  The  results  in  Figure  4 are  interesting  since 
one  would  expect  a decrease  in  frequency  as  damping  is  increased.  Stiffen- 
ing effects  due  to  damping  are  found  elsewhere  [5]  and  could  be  due  to  the 
manner  in  which  it  is  introduced  in  the  problem.  Figure  4 also  shows  that 
the  results  for  the  fully  damped  beam  are  approached  rapidly  as  damping 
is  increased.  Figure  5 shows  the  vibrating  frequencies  on  a complex 
plane  for  a beam  with  rotation  flexibility  and  damping.  The  load  is  zero 
(non-rotating  beam)  and  the  rotation  spring  is  kept  at  k2  = .1  while  the 
dashpot  value  changes.  The  arrows  adjacent  to  the  curves  show  the  direc- 
tion of  the  values  on  the  curve  as  damping  increases.  For  zero  damping, 
the  results  are  purely  imaginary  and  are  approximately  .54  and  15.5  for  the 
1st  and  2nd  frequencies.  As  damping  increases,  the  frequencies  become 
complex  with  the  imaginary  components  decreasing  for  the  first  eigenvalue 
and  increasing  for  the  second.  The  behavior  of  the  first  frequency  is 
interesting.  As  the  damping  value  increases,  the  imaginary  component 
vanishes  (as  also  seen  in  Figure  4)  as  if  the  system  becomes  critically 
damped.  However,  the  real  component  can  also  be  followed  on  the  complex 
plot  and  the  beam  appears  to  vibrate  again  in  this  first  mode  as  damping 
increases  further.  For  sufficiently  large  damping  the  frequencies  seem 
to  approach  those  for  a beam  which  is  damped  at  the  support.  Points  on 
the  real  axis  represent  zero  motion  but  move  with  changes  in  damping 
values  to  points  on  the  real  axis  where  it  is  intersected  by  a branch  or 
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mode.  Figure  6 shows  the  same  results  for  load  Q = 25.  A final  result  for 
rotation  flexibility  is  shown  in  Figure  7 for  k2  = Here  the  rotation 

spring  is  relatively  stiff  and  the  fully  damped  results  are  rapidly 
approached  with  initial  effects  for  near  zero  dashpot  values  overshadowed. 

The  effect  of  support  damping  on  the  frequencies  on  beams  with 
deflection  flexibility  are  shown  in  Figures  8-10.  A decrease  in  frequency 
with  increased  damping  is  seen  here.  Figures  9 and  10  show  the  effect  on 
the  complex  plane  for  Q = 0 and  Q = 200,  respectively. 

Finally,  Figure  11  shows  the  results  of  a case  when  both  rotation 
spring  and  deflection  spring  flexibilities  are  allowed.  Little  effect 
on  frequency  is  noted  as  the  rotation  spring  is  varied  while  the  deflec- 
tion spring  and  dashpot  remain  constant  in  value.  An  almost  parallel 
increase  in  frequencies  are  found  when  the  rotation  dashpot  is  increased 
in  value.  The  investigation  into  the  response  of  the  beam  with  all  springs 
and  dashpots  finite  was  limited  to  those  shown.  A study  should  be  per- 
formed to  indicate  areas  where  the  effects  on  beam  response  will  be  most 
pronounced. 
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Figure  9.  Effect  of  Support  Damping  on  Frequencies 
of  Beams  with  Deflection  Flexibility:  Real  vs 
Imaginary  Part  of  Eigenvalue 
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Figure  11.  Effect  of  Rotation  Spring  on 
Frequencies  for  Fixed  Finite  Deflection 
Spring 


Q = 25.0 
kx  = 25.0 
Cj  = 5.0 
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